Exploration of resonances by analytic continuation in the coupling constant 
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The energy and the width of resonance states are deter- 
mined by analytic continuation of bound-state energies as a 
function of the coupling constant (potential strength). The 
advantage of the method is that the existing techniques for 
calculation of bound states can be applied, without any mod- 
ifications, to determine the position of resonances. Various 
numerical examples show the applicability of the method for 
three-body systems, including the excited states of the 6 He 
and 6 Li. 

PACS numbers:25.70.Ef, 21.45. +v, 27.20. +v 



Although many of excited states are resonances JjJ] 
in nuclear physics, they seldom receive special attention 
in theoretical nuclear structure calculations. The rea- 
son is obvious: The solution of the nuclear many-body 
Schrodinger equation poses insurmountable difficulties 
even for bound states in most of the cases, and the so- 
lution for the resonances is mostly limited to two-body 
problems. The theoretical nuclear model calculations 
therefore usually treat the resonances as bound states. 
This approximation, however, can only be justified for 
very narrow resonances. 

In the recent years many theoretical and experimental 
nuclear physicists have concentrated their research ac- 
tivity on unstable nuclei and the new phenomena discov- 
ered have revived the interest in resonance states as well. 
These nuclei are very rich in resonance states, and one 
of the most spectacular phenomena, the "Borromean" 
binding 0j , is also closely related to resonances. In a Bor- 
romean system (e.g., 6 He or n Li) the bound three-body 
system has no bound binary subsystem, but the two-body 
subsystems usually have a resonance. Besides these two- 
body resonances, one can find an example for three-body 
resonances and for even more complicated ones as well 



Although there are elegant and powerful techniques to 
find resonances as poles of the S-matrix (see e.g. jL|), 
these are mostly limited to simple interactions and sys- 
tems. For composite nuclear systems the theoretical ap- 
paratus has been formulated for finding the energies of 
bound states. The easiest way is thus to use the machin- 
ery developed for bound states to gain information on 
resonances. There are indeed two methods widely used 
in nuclear and atomic physics which rely on bound state 
type solutions, that is, use square integrable functions. 
The first, the "real stabilization method" solves the 
Schrodinger equation in a box and exploits the fact that 



the energy of the resonance state remains stable against 
the change of box size within certain limits (the wave 
function is localized inside the box), while the energies 
of the continuum states rapidly change. The second, the 
complex scaling method (CSM) 0, rotates the coordi- 
nate r by e l6 r and transforms the continuum resonance 
wave function to the normalizable wave function of the 
bound state. Both of these methods have the advantage 
that they use square integrable functions and therefore, 
after some modifications, the already existing methods of 
bound-state type can be applied. These two methods are 
known to be in an intimate connection [§J . A slight disad- 
vantage of the real stabilization method is that in order 
to determine the width of a resonances one has to solve 
the Schrodinger equation with many different box sizes 
and that it becomes computer time consuming for com- 
posite systems. Moreover, the box for a many-particle 
system is a many-dimensional object and sometimes it is 
difficult to find the appropriate intervals to be changed. 
In practical cases the method turns out to be suitable 
for narrow resonances only. The complex scaling method 
has proved to be capable of exploring wider resonances. 
Its computational burden is the calculation of complex 
matrix elements and the solution of complex eigenvalue 
problem. The source of the real hardship is that the 
Hamiltonian is nonhermitian and the variational princi- 
ple to find an energy minimum cannot be applied. A 
useful guiding principle for selection of the basis, there- 
fore, has been lost and it becomes complicated to choose 
the appropriate trial functions. 

Kukulin et. al. suggested an even simpler method to 
study resonant states; the analytical continuation in the 
coupling constant (ACCC). Their approach is based on 
the intuition that the resonances can be thought as the 
continuation of bound states when the attraction of the 
interaction decreases. More precisely, they analytically 
continue the energy of the bound state as a function of 
the strength of the potential to the complex plane to re- 
veal the width and energy of the resonance. The authors 
parameterize the square root of the energy by a Pade ap- 
proximation as a function of the coupling constant of 
the potential (potential strength). The coefficients of 
the Pade approximation are determined by solving the 
Schrodinger equation for the coupling constants giving 
bound states. The computational demand thus amounts 
to solving the Schrodinger equation for bound states for 
several different potential strengths. The Pade approxi- 
mation is chosen to approximate the square root of the 



1 



energy as a function of the coupling constant because it is 
more general and more powerful than other possibilities. 
For example, unlike a Taylor expansion it can simulate 
singularities of a function near the threshold. 

The ACCC method was hitherto applied in a few sim- 
ple test cases for two-body resonances e.g., the analyt- 
ically solvable square-well potential and a macroscopic 
a-a system with a simple potential. The application was 
limited because to fix the coefficients of the Pade ap- 
proximation one has to solve the bound state problem 
accurately. Nowadays, there are several reliable meth- 
ods at hand to solve few-body problems. Before going to 
some really challenging applications, it is unavoidable to 
test the performance and to learn the limitations of the 
ACCC method. 

The aim of this paper is to show that the ACCC is re- 
ally a powerful method to study the resonances of nuclear 
systems. A mathematical proof is not available, and the 
analytical solution is available only for a few very special 
cases. Therefore we must rely on numerical examples 
for various cases. We compare the results of ACCC to 
other methods such as direct numerical integration(DNI) 
or CSM. 

The ACCC assumes that the Hamiltonian of the sys- 
tem is written as H(X) = Hi + XH2 , where H2 is the 
attractive part of the interaction. By decreasing A the 
bound state approaches the threshold (at Ao the energy 
is -E(Ao) = 0) and may become a resonance or a vir- 
tual state. In Ref. Q it has been shown that for a 
two-body system, near the threshold the square root of 
the energy behaves as ki (A) ~ y/X — Xq for I > and 
feo(A) ~ (A — Ao) for S'-wave. Introducing a variable 
x = y/X — Ao the analytic function ki has two branches 
ki(x) and ki(—x). Due to the analyticity of these func- 
tions we can continue them into the resonance region 
(A < Ao) from the bound states (A > A ). Motivated by 
the above functional form of ki near the threshold, the 
Pade approximation of the form H 
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is used for the analytical continuation. In principle Co is 
to be zero but reserved to take care of possible errors in 
the determination of Ao- In practice Co is found to be 
much smaller than other c values. 

The coefficients of the polynomials are calculated in 
the bound state region and therefore are real. If A < Ao 
x is pure imaginary and ki(x) becomes complex. The 
energy and width of the resonance state is given by 
E — iT/2 — ki 2 To determine the coefficients we are to 
solve the bound state problem for various coupling con- 
stants A (> Ao) and try to find the threshold value Ao 
(fcj(Ao) =0). To have a reliable approximation one has 
to know the accurate values of the coefficients in the Pade 
approximation, that is, one has to solve the bound state 
problem to high (typically 4 or more digits) accuracy, 
especially at the threshold. The stochastic variational 



approach |lO|]ll| ] seems to be quite appropriate for this 
purpose. 

Since one has to solve the bound state problem many 
times to determine the coefficients, one may think that 
the ACCC would be computer time consuming for larger 
systems. Due to the simple linear dependence of the 
Hamiltonian on the coupling constant, however, one does 
not have to recalculate the matrix elements and the com- 
putational load is just the rediagonalization of the Hamil- 
tonian for different values of the coupling constant. This 
property makes the application especially simple. 

To solve the bound state problem we use the stochastic 
variational method. Of course any other method may be 
suitable. A point to be emphasized here is that by chang- 
ing the coupling constant one goes from deeply bound to 
weakly bound states and finally to the threshold. All 
these states should be accurately treated and in a varia- 
tional calculation, for example, one has to choose a basis 
which adequately spans the configuration space. 

To illustrate the method we start with a simple exam- 
ple. We consider a system of two-particles with mass m 
interacting via a two-range Gaussian potential: 



V(r) = -8A exp[-(r/2.5) 2 ] + 2 exp[-(r/5) 2 



(2) 



where Ti = c = m = 1. This simple problem can be eas- 
ily solved by DNI [113]. The trajectories of resonances of 
different partial waves by ACCC and by DNI are com- 
pared in Fig. 1. For simplicity we assumed M = N 
in Eq. (]l|). To determine the 2M + 1 coefficients one 
has to solve the bound state for 2M + 1 different cou- 
pling constants (Ai, . . . , A 2 a/+i). Once the bound state 
energies are known, the coefficients of the Pade approx- 
imation can be extracted by solving a system of linear 
equations. Due care must be taken to avoid numerical 
problems. The ACCC results agree very nicely with the 
DNI results including the IS excited resonance state. 

The second example tests the effectiveness of the 
method for wide resonances. The wide resonances may 
cause serious difficulties in many cases. In this example 
we have used the same potential as Eq.(pl) but with a 
lower barrier (Vb = 0.25 instead of Vq = 2). The results 
of ACCC and DNI agree very well (see Fig. 2). The re- 
sults are surprisingly good considering the fact that, in 
the last few points where the attraction is very weak, the 
widths of the states are two times of their energies (one 
may not really call them resonances). 

The third example is a three-body case. Three bosons 
of mass, h 2 /m — 41.47, interact via the potential 



V(r) = -120 exp[-r 2 ] + 3A exp[-(r/3) 2 



(3) 



The energy and length are in units of MeV and fm. Note 
that there is no two-body bound state. We compare the 
results of ACCC to those of CSM. This simple example 
is selected because the CSM might have some inaccu- 
racy in more complex cases, and we want to address the 
applicability of ACCC in a clean test case. The three- 
body bound state problem for zero total angular mo- 
mentum has been solved by the stochastic variational 
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method by using a Gaussian basis O). The trajecto- 
ries of ACCC and CSM as a function of A are shown in 
Fig. 3. The two methods give the same resonance posi- 
tion for a wide range of potential strength, justifying the 
suitability of ACCC. Pade approximations with M=3, 5 
and 7 are used. By increasing the number of terms the 
agreement slightly improves but already M=3 terms give 
good results. Except for the threshold value, one has the 
freedom to choose those coupling constants for which the 
bound state problem is solved to determine the coeffi- 
cients of the Pade approximations. The position of the 
resonance, of course, depends on the choice of the set of 
A. We found that if the values of coupling constants used 
are distributed into a wide range the dependence of the 
resonance parameters E and T on the input values are 
relatively small. The sensitivity can be easily controlled 
by comparing the results starting from several adequately 
chosen input sets. 

The resonance states of the a + N + N three-body 
model are chosen as a practical example. These reso- 
nances have been very intensively studied in the past 
and the investigation is continued with great elan 
Previously, we described the bound states of the 
^le and 6 Li nuclei in an a + n + n and a + p + n-type 
microscopic cluster model with the stochastic variational 
method flllliul . To determine the resonances of these sys- 
tems by ACCC, we used the same type of bound state 
solutions. The model assumes an alpha cluster but uses 
a fully microscopic six-body wave function. The nuclc- 
ons interact via the Minessota effective interaction (sum 
of the Coulomb and spin-isospin dependent central and 
spin-orbit potentials) |L5| . The strength of the attraction 
of the potential is controlled through the space exchange 
mixture parameter u. Note, that due to the simple har- 
monic oscillator shell model description of the alpha par- 
ticle, the change of the potential parameter u does not 
change the energy of the a-particle and thus the three- 
body threshold remains the same. The wave function 
of the system is taken as linear combination of terms 
describing the a(NN) and (otN)N arr ang ements. The 
details of the model is given elsewhere [[L4j . To compare 
the results to other calculations the parameters and the 
model space of Ref. jj| have been adopted. 

Table I compares the results of ACCC with those of 
CSM. The CSM is not an exact solution, and especially 
for complex systems, the results of CSM might have in- 
herited some inaccuracies from the underlying gaussian 
expansion and numerical inaccuracies due to the complex 
arithmetics ||. The results of ACCC and CSM therefore 
should be consistent within "error bars" . Bearing in mind 
that both methods attempt to solve a composite system 
starting from a fully microscopic model, the agreement 
of the result can be considered to be very good. These 
examples show that the ACCC can be combined with a 
microscopic structure model. The ACCC might be ap- 
plied in combination with other microscopic method like 
shell model or Hartree-Fock method. 

In a three-body system, the form of singularity in the 



coupling constant is likely to be different from that of the 
two-body case. The need for higher order terms in the 
Pade-approximation may reflect this incorrectness. The 
knowledge of the analytical form of the near threshold 
singularities would make the convergence (in the terms 
of Pade-expansion) faster. 

In summary, we have shown through various exam- 
ples that the ACCC is really a powerful method to cope 
with resonances of nuclear systems. The application of 
the ACCC method has been made possible by the re- 
cent developments in solving bound state type problems 
and the increase of computational power. The results of 
the ACCC method have been compared to those of other 
solutions and are found to be in good agreement. 

The advantage of the ACCC over the CSM is that 
one does not need to calculate complex matrix elements, 
that one does not have to use complex arithmetics on 
the computer. Unlike the real stabilization method, one 
does not need to recalculate the Hamiltonian on differ- 
ent bases many times. On the other hand, to have a 
reliable solution by the ACCC method, one has to solve 
the bound state problems very accurately. This require- 
ment may pose certain limitations in applications. The 
determination of resonances by the direct solution of the 
Faddeev-equation |ll||l7j is certainly superior, but the 
ACCC seems to be more easily applicable in the frame- 
work of microscopic models. 

The results encourage the application of ACCC to 
study the resonance states of few-nucleon and few-cluster 
systems. 
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FIG. 1. The resonances for the two-range Gaussian poten- 
tial of Eq.(2). The ACCC results with M=iV=5 are com- 
pared with those of direct numerical integration (DNI). The 
trajectories of the resonances are plotted as a function of the 
coupling constant A at intervals of 8X = 0.04. 



FIG. 2. The trajectories of wide 5-wave resonances. The 
same as Fig.l but the barrier height is reduced as described in 
text. The Pade approximation with M—N=19 is used. The 
energy and the width of the resonance are also shown as a 
function of A. 
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FIG. 3. The resonances of the three-boson system inter- 
acting via the potential of Eq. (3). The ACCC results with 
M=N=7 are compared with those of the complex scaling 
method (CSM). The trajectories of the resonances are plot- 
ted as a function of the coupling constant A at intervals of 
8X = 0.1. 



TABLE I. Comparison of resonance energy and width be- 
tween ACCC and CSM. The energy is from the three-body 
threshold. 
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